Bell test with time-delayed two-particle correlations 
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Adopting the frame of mesoscopic physics, we describe a Bell type experiment involving time- 
delayed two-particle correlation measurements. The indistinguishability of quantum particles results 
in a specific interference between different trajectories; the non-locality in the time-delayed corre- 
lations manifests itself in the violation of a Bell inequality, with the degree of violation related to 
the accuracy of the measurement. In addition, we demonstrate how the interrelation between the 
orbital- and the spin-exchange symmetry can by exploited to infer knowledge on spin entanglement 
from a measurement of orbital entanglement. 
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I. INTRODUCTION 

Fundamental quantum phenomena, such as non- 
locality and entanglement of quantum degrees of free- 
dom, have regained a lot of interest recently, mainly due 
to their potential usefulness as a computational resource. 
Mesoscopic physics provides a new platform for the in- 
vestigation of these phenomena, important issues be- 
ing the creation, quantification, and verification of non- 
locality/entanglemcnt. In this paper, we describe an ex- 
periment where two electrons with different orbital wave 
functions are superposed in an interferometer and ana- 
lyzed in a Bell type experiment involving two-particle 
correlation measurements, see Fig. [1] The particular 
feature of this Bell test is the replacement of the four 
different settings of local detectors in the original setup 
by four different time-delays in the measured correlators. 
The main physical property we want to exploit is the in- 
distinguishability of quantum particles, which results in 
a specific interference between different trajectories. We 
wish to convey three messages: first, the non-locality in 
the time-delayed correlations due to indistinguishability 
manifests itself in the violation of a Bell inequality. Sec- 
ond, the degree of violation is related to the accuracy 
of the measurement and is reduced, once the local mea- 
surement can distinguish between the different orbital 
wave functions of the particles. The above two items 
refer to spinless objects (or particles with equal spin). 
Third, adding the spin degree of freedom, we show how 
the symmetry relation between spin- and orbital compo- 
nents allows to extract information on spin-entanglement 
from an orbital measurement. 

By now, numerous proposals have been made how to 
create entangled states in mesoscopic setups, both for 
orbital- and spin degrees of freedom^. The verification 
and quantification of entanglement can be carried out 
using Bell inequality checks 2 or state tomographj*2ii£. 
The indistinguishability of (spinless) particles produc- 
ing a two-particle Aharonov-Bohm effect and entangle- 
ment has been exploited in a Hanbury-Brown Twiss 
interferometer^; here, we offer an alternative imple- 
mentation which makes use of an electron beam split- 
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FIG. 1: Particles incident in leads u and d with wave functions 
f(x) and g(x) are mixed in a four-terminal splitter (character- 
ized by the mixing angle 9) and analyzed at xi = X2 through 
measurement of time-correlations during the time intervals 
ti £ A and (266 (time axis drawn perspectively into the 
plane). We are interested in the entanglement of the lead 
indices u and d with respect to bipartitioning of the system 
between the time intervals A and B. For details on the im- 
plementation of the source (shaded area) see Fig. [4] 



ter. Using the same setup as discussed here, Burkhard 
et a/4 have demonstrated how to distinguish between 
singlet and triplet spin-states by a measurement of zero- 
frequency cross-correlations. Below, we exploit that the 
orbital measurement of the Bell parameter preserves the 
spin-entanglement; this feature allows us to find a lower 
bound on the concurrence of the spin wave function. 

In the following, we consider a setup with two incom- 
ing leads, denoted as u and d, connected to two outgoing 
leads u and d through a reflectionless four-terminal beam 
splitter, see Fig.[T] At time t = two electrons with nor- 
malized orbital wave functions f(x) and g(x) and com- 
mon spin state x(°i> ^2) are injected into the leads u and 
d. The state, factorizable in orbital and spin parts (with 
factorized orbital part and general spin part), is conve- 
niently written within a second quantized formalism, 



|*in) = J dxidx 2 f(x 1 )g(x 2 ) 



(1) 
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here, ^ aa (x) creates electrons at the position x in lead 
a and |0) is the vacuum state with no electrons. After 
mixing in a four-terminal splitter (with mixing angle 9), 
we will analyze correlations in the system through detec- 
tion of particles in time (see Fig. [1]) or space separated 
intervals A and B (see Fig. [3]). Our focus then is on the 
entanglement of the lead indices u and d with respect to 
bipartitioning of the system between the time or space 
intervals A and B. 



II. BELL TEST 

The particles in the outgoing leads u and d are sub- 
jected to a Bell test expressed through time-resolved 
current-current correlators in the leads a\ and a.2, 
ai, a<x € {u, d} (both auto- a% = a 2 and crossed- ai ^ a 2 
correlators are considered), 



dti / dt 2 (I ai (xi,t 1 )I a2 (x 2 ,t 2 )}, 



C ai a 2 (AB) — — 

Al ~ J i ./ ;; 

(2) 

where I a (x, t) is the total current operator (summed over 
spin degrees of freedom) in lead a at position x and time 
t. The time integration is taken over a finite time interval 
A=[t A - St/2, t A + St/2] (same for B) with the width 
St accounting for the finite time-resolution of the current 
measurement, cf. Fig.Q] the limit <5£ — > corresponds to 
a measurement of the instantaneous current. In the fol- 
lowing, we will assume that all correlators are measured 
at some fixed symmetric position x\ = x 2 and omit the 
coordinate variable. 

With only two electrons present in the system and for 
non-overlapping time-intervals AnB = 0, the correlation 
function C aiCC2 (AB) is proportional to the joint probabil- 
ity Pa ia2 {AB) for the detection of two particles during 
the time intervals A and B in the leads a,\ and a 2 , see Ref. 
9. There are four distinct possibilities to distribute two 
electrons between the outgoing leads and we can define 
the properly normalized E a a P aiol2 (AB) = 1) proba- 
bilities as 



Bcx±cx2 

(AB) 



C'a 1 a 2 ( J 4J3) 

(AB)' 



(3) 



Out of these, we define the two-particle Bell inequality 
in the Clauser-Horne^ form in the same way as it is 
done in the usual optics context^: we introduce the Bell 
correlation functions 

E AB = [Pun - Pud - Pdu + Pdd]AB (4) 

and obtain the Bell inequality 

\E A b -E A b> +E a >b +E A 'B'\ < 2. (5) 

Here, the polarizations ± in the optics context are re- 
placed by the lead indices u and d and the role of the four 



different polarization settings of the detectors is played 
by four different time intervals A, B and A', B' . The 
violation of this inequality for a particular choice of time 
intervals shows that non-local correlations are present in 
the system, i.e., the result of the measurement cannot be 
simulated by any local-variable theory. 

Let us demonstrate that the above Bell inequality in- 
deed can be violated by the incoming state |T]) after 
proper projection. We then have to calculate the four 
current-current auto- and cross-correlators C Q1Q2 with 
ai = a 2 and a\ ^ a 2 , respectively. This is done within 
the scattering matrix approach to quantum noised: we 
assume that the Fourier components f(k) and g(k) of the 
single-particle wave functions are concentrated near the 
wave vector ko > 0, allowing us to linearize the energy- 
momentum dispersion near fco- The time evolution of the 
incoming state ([T]) then is described by the propagation 
of the single-particle wave packets f(x) and g(x) with 
constant velocity vq = hko/m to the right, f(x, t) = /(£) 
and g(x,t) = <?(£), where £ = x — v$t is a retarded vari- 
able. With the scattering matrix of the beam splitter 
(parametrized by the angle 8), 



cos f — sin t 
sin 6 cos 9 



(6) 



we can express the current operators I a (x,t) in the out- 
going leads a 6 {u, d} through the electronic scattering 
states. Averaging the product of current operators in 
Eq. @ over the incoming state l^in) one arrives at the 
results 

(L(^)iu(^)) = (ev ) 2 {(co8 2 8\f(^)\ 2 (7) 
+ sm 2 e\g(&)\ 2 )S(Z 1 -&) 
+ sin 2 0cos 2 0[|/(e 1 )| 2 | 5 (6)| 2 + |.g(6)| 2 |/(6)| 2 

-Q(/(ei)9*(6)5(6)r(6)+c. c .)]}, 

(4(£iKd(6)) = (e«o) 2 {cos 4 e|/(ei)| 2 |5(e 2 )| 2 (8) 
+ sin 4 %(£ 1 )| 2 |/(6)| 2 
+ sin 2 cos 2 0g (/(£i)g*(£i) 5 (6).r(6) + c.c.)}, 

where Q = J2 ail72 x( cr i, cr 2 )x*(cr 2: <J\) describes the over- 
lap between the spin states of the two electrons in the 
incoming state ([1]). The two other correlation func- 
tions (Jd(£i)/d(6)> and (£l(£i)iu(6)) are obtained by 
exchanging cos 9 and sin# in Eqs. ([7]) and (jHJ). Substi- 
tuting these expressions for the current correlators into 
Eq. @ and integrating over (non-overlapping) time inter- 
vals A and B, one arrives at the Bell correlation function 
Eq. Q 

E AB = - cos 2 (2fl) - sin 2 (2g) Q ^ + g^L , (9) 

where we have introduced the particle densities Fa,b and 
Ga,b averaged over the time intervals A and B, 



F ^ = Tt 



dt\m)\ 2 



(10) 



t£A,B 



(and similarly for Ga.b with / replaced by g). The over- 
lap Sa,b between different single particle wave functions 
reads 



Sa,b = y + I #/(0s*(0- 
01 JteA.B 



(11) 



In the following, we apply the result ([9]) first to spinless 
fermions and plane wave states / and g and confirm the 
violation of the Bell inequality in this simple situation. 
We then proceed with a rederivation of the expression © 
with space-like separated measurement intervals A and 
B in order to make the origin of the entanglement more 
transparent. A formulation in terms of reduced density 
matrices leading to an expression of the Bell correlator 
in terms of concurrences completes the discussion. 



A. Spinless fermions 

We first concentrate on spinless particles; this situa- 
tion can be realized by preparing the two electrons in 
equal spin-states, xtVi) ^) = ^o-iT^T with correspond- 
ing overlap Q = 1. To begin with, we choose a plane 
wave form for the wave packets with different momenta 
k\ and ki close to fco ( m order to allow for the linearized 
spectrum), f(x) — exp(ifcia;) and g{x) — exp(ifc2x). The 
correlation function Eab takes the form 



E A b = -cos 2 (26>) - V sin 2 {28) cos^ AB , 



(12) 



where lpab — Soj {JLa ~ is) is the relative phase shift 
accumulated by the two waves between the two measure- 
ment intervals and Soj = Vo(ki — kz) is the frequency 
mismatch between the two plane waves. The phase shift 
ifiAB replaces the angle between the two polarizers in the 
conventional Bell setup. The visibility factor < V < 1 
accounts for the width St of the time interval, 



V = 



sin 2 (SujSt/2) 
{5ujSt/2f ' 



(13) 



The other correlation functions involving intervals A' and 
B' are obtained in the same way; their combination into 
the Bell inequality Eq. ([5]) produces a maximal viola- 
tion for the angles (Pab = Va'b = PA'B' = 7r /4 and 
y>AB' — 37r/4, corresponding to measurement intervals 
with relative distance ts = tA + t/8, Ia 1 = tA+T/4, and 
tb' = tA + 3t/8, where r = 2it j ' Suj and tA is an arbitrary 
reference time. In this case the Bell inequality Eq. (J5]) 
reduces to 



£ = \vV2sm 2 (28) + cos 2 (20)| < 1. 



(14) 



For given V, the maximal violation £ max = V \[2 is 
reached for a symmetric beam splitter with 8 = tt/A, 
cf. Fig. [2] Furthermore, the maximally allowed degree 
of violation £ = \[2 can be attained only for V ~ 1, 
corresponding to a short time measurement of the cur- 
rent value with St < 1/voSk — Ti/Se: hence, the maximal 
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FIG. 2: Bell inequality violation for maximal visibility V — 1 
versus mixing angle 9 (solid line: spin-triplet states; dashed 
line: spin-singlet state). 



violation of the Bell inequality can be obtained for indis- 
tinguishable particles, while a time interval with length 
beyond Heisenberg's uncertainty bound St > fr/Se allows 
for a distinction between the two particles and the Bell in- 
equality cannot be violated in this classical situation. For 
V < 1/ V% the Bell inequality is always satisfied. Within 
the region 1/V2 < V < 1, the Bell inequality ([H]) is 
always violated for any mixing angle < 8 < tt/2, al- 
though to a lesser degree then in the symmetric point 
8 = tt/4. 

Let us discuss the physical origin of the violation. The 
correlator Eq. © measured in the Bell test is finite, pro- 
vided that both electrons are detected within the time 
windows A and B; in this case, it is proportional to 
the probability P ai0t2 (AB). Although, formally, the elec- 
trons have different energies e and thus are distinguish- 
able in principle, given a small time resolution St < h/Se 
of the local current measurements one cannot distinguish 
between the energies e\ = Tiv^ki and £2 = ftvofe- Under 
this circumstances the electrons indeed can be consid- 
ered as indistinguishable particles. Then, according to 
the rules of quantum mechanics, there are two quantum 
alternatives contributing to a coincident detection of the 
electrons in A and B: either the electrons with energies £\ 
and £1 are detected in the time windows A and B, respec- 
tively, or vice a versa. These two alternatives contribute 
to the measurement outcome with different phases: in the 
first case, the phase factor acquired by the two-particle 
wave function after the first measurement at tA due to 
the propagation of the second particle until ts is given by 
exp[— i£2{tB~tA)], while in the second case this phase as- 
sumes the value exp[— i£\ (ts —tA)]- The phase difference 
between the two alternatives leads to quantum interfer- 
ence and a corresponding oscillatory dependence (with 
frequency Suj — (£2—£i)/fo) of the probability P ai a 2 {AB) 
as a function of time, with an amplitude proportional 
to the visibility factor V. The precise bound on St al- 
lowing for a violation of the Bell inequality is given by 
1/V2 < V < 1 or (e 2 — £i)St < 2K, corresponding to a 
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measurement where the Heisenberg uncertainty principle 
for energy-time variables is violated. 



B. Space-separated domains 

In order to understand better the nature of the entan- 
glement observed in (114|) . we consider a slightly different 
experiment, where instead of using time-separated de- 
tection intervals, the two observers Alice and Bob are 
measuring the simultaneous appearance of particles in 
spatially separated regions A and B of the setup, see Fig. 
[3J for particles with a linear dispersion, these two exper- 
iments are equivalent since a time delayed measurement 
with St = t 2 — t\ at the point x corresponds to a coinci- 
dent measurement at time t with Sx — (t 2 — tijVQ. We 
first concentrate on plane-wave incoming states, where 
the present setup with spatially separated detectors pro- 
vides additional insights. In particular, we will see that 
it is the projection of the non-entangled incoming state 
onto the two domains A and B that defines a biparti- 
tion of the system with respect to which the lead index 
becomes entangled. On the other hand, in order to per- 
form a Bell inequality check, we need a set of local 'rota- 
tions' of the measurement apparatus: in our setup, the 
parameters generating a suitable set of local 'rotations' 
are determined by the distance between the measurement 
domains A and B and by the mixing angle 6. 

A central element in our discussion below is the in- 
terchangeability of mixing U ® U and projection Vab 
onto the domains A and B, where U denotes the one- 
particle scattering matrix of the beam splitter and the 
tensor product U®U acts on our two-particle state. This 
interchangeability is a trivial consequence of these two 
operations affecting different degrees of freedom, coordi- 
nates X\ and X2 and lead indices u and d. In terms of 
these operators, we can relate the incoming and outgoing 
states via 



in u i) 



VabU®UW 



(15) 



Assuming that the projections onto A, B in the outgoing 
leads and onto A, B in the incoming leads (see Fig. [3j 
are ballistically separated (i.e., the measurement in A, 
B involves the appropriate ballistic delay time) we can 
write 



\*ab) 



U(S)UVar\^ 



(16) 



Hence, in our discussion we are free to interchange the 
two operations of mixing and projection. 

Consider then an incoming state (before mixing) with 
single-particle wave functions f(x) = e zklX and g(x) = 
e zk 2 x w jt n shifted momenta. The state incident from 
leads u and d can be written as a simple Slater deter- 
minant, 
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FIG. 3: Particles incident in leads u and d with wave functions 
f(x) and g(x) are mixed in a four-terminal splitter and ana- 
lyzed through measurement of equal-time correlations within 
the space intervals A and B centered around xa and x b ■ The 
interchangeability of projection (to the intervals A and B) 
and mixing allows to shift the measurement intervals to the 
positions A and B in the incoming leads; provided that the 
measurements in A and B and in A and B are ballistically 
delayed in time, the measurement outcome is the same. 



and thus is non-entangled. The lead index x £ {u, d} of 
the electron field operator ip x is conveniently regarded as 
a pseudo-spin. 

To start with, we analyze the coincident detection of 
two particles within the non-overlapping regions A and 
B of the incoming leads, see Fig. [3J and select only 
those events, where each of the observers (Alice in A 
and Bob in B) finds only one particle. For two parti- 
cles, this projection can be described by the operator 
V Ab — N(A)N(B), with the particle number operator 

N(X) = j x dx (ipt(x)i{>u(%) +i>\(x)4>d( x )) counting par- 
ticles in the region X of the incoming leads. Projecting 
the incoming state |T7J) one arrives at the state 



I* 



AB, 



= J dxidx 2 [fA( x i)9s( x 2) 
+f s ( Xl )g A (x 2 )}$i(xSi(x 2 )\0), 



(18) 



* m ) = / dx 1 d X 2f(x 1 )g(x 2 ) 1>l(xi)mx 2 )\0), (17) 



where fx{x) and gx{x) are equal to f(x) and g(x) for 
x G X and vanishing outside. This projected state is no 
longer a simple Slater determinant and describes a two- 
particle state entangled in the lead indices and shared 
between the regions A and B of the incoming leads. It 
is instructive to rewrite the state (|18p in a pseudo-spin 
notation: Assuming for simplicity that the intervals A 
and B are reduced to individual points Xa and x§ we 
have 

\*%) oc e iVAs/2 \]) A \[) s + e-^^ 2 \l) A \V s , (19) 

where | |)jc an( l I i)x denote states of particles localized 
in X and residing in lead u and d, respectively; the orbital 
part of the wave function contributes the phase factors 
exp(±i(^ j 4^/2) with (fiAB = ^{ x a~ x b)^ where Sk = k\ — 
k 2 is a momentum mismatch. The projected state (|19l) is 
in fact maximally entangled in the lead- or pseudo-spin 
index with respect to bipartitioning the system between 
the regions A and B. In the following, we wish to detect 
this entanglement in a Bell test. 
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The implementation of a Bell test relies on the ability 
to locally change the pseudo-spin basis of the particles. 
To do so, we transmit the original incoming state Eq. 
(jTTJ) through a beam splitter before measuring the pres- 
ence of particles in the intervals A and B, now located in 
the leads u and d to the right of the mixer, see Fig.O the 
mixing then acts as an equal rotation of the (pseudo-spin) 
basis u, d for both particles. However, such a global rota- 
tion of the original basis alone is not sufficient to perform 
the Bell test, as locally distinct rotations are required as 
well; the latter are implemented through different choices 
in the separation Sx = x 2 — %\ between the regions A and 
B. Exploiting the interchangeability of projection and 
mixing, cf. Eqs. (|15|) and (fT6|) . we see that this change 
in distance results in a relative rotation with the angle 
'PAB = <Pab around the original (u, depolarization axis 
of the pseudo spins, see Eq. (fl9|) . Writing the outgoing 
state (|16[) in pseudo-spin notation, we obtain the expres- 
sion 

= -cos(^W2)sin(2^ lT)s ~ llUli)B 



+ cos((^4b/2) cos(26>) 
+ism(if AB /2) 



V2 

1)a\1)b + \1)a\1)b 

T)aU)b-U)a|T)b 



V2 



(20) 



This projected state describes two spatially separated 
localized particles with entangled pseudo spin indices. 
Choosing different space separations between the regions 
A and B allows one to change the phase pab and mix- 
ing by U <£> U generates a second rotation parametrized 
by the angle 6. Calculating the joint probabilities 
Paia 2 (AB) oc \(aia 2 \^ A u B )\ 2 for the four settings a\a 2 € 
{uu, ud, du, dd} we find the Bell correlation functions 
Eab as given by ^ and choosing appropriate angles 
Vab, *Pab', fA'B, fA'B> an d 6 one finds the Bell inequal- 
ities violated. 



C. Density matrix formulation 

In a last step, we reformulate our analysis in terms of 
density matrices and express the Bell inequality in terms 
of concurrences of density matrices reduced after projec- 
tion to the intervals A and B. We rewrite the projected 
state (TT5)) incident from leads u and d in pseudo spin 
representation, |\l/x§) = / dxxdx 2 Y& AB {x\,x 2 )), wnere 



\^ab(xi,x 2 )) = f A (xi)g s (x2) | T)aII)j 
+/s(^i)5a( 2; 2)|I)^|T)S' 



(21) 



The joint measurement of the pseudo-spin index in the 
regions A and B is described by the coordinate-reduced 
two-particle density operator 

Pab « / dxidx 2 \$! Ab{ x Ii x 2)){^ Ab{ x Ii x 2)\- 



We introduce the two-particle pseudo-spin basis {|TT)j 
|TI); lit); III)}! where the first (second) arrow refers to 
the particle localized in A (B); the normalized density 
matrix then assumes the form 



Pab- 



F A Gt 



-Ga f b 



/0 







f aGb 
~S\Sb 





—Sa^b 
Ga f b 








where F x — J x \f(x)\ 2 dx, Gx — J x \g{x)\ 2 dx and Sx — 
f x f(x)g*(x)dx with X £ {A,B}. Although initially 
the two particles have been in a pure state, the reduced 
density matrix (f2"2")) corresponds to a mixed state with 
p\b ^ Pab- The calculation of the entanglement in the 
mixed two-particle state Eq. (|22[) corresponds to finding 
the concurrence C(pab) °f a two-qubit problem, and thus 
can be calculated following the scheme introduced by 
WoottersH, C(p^ B ) = max{0, V\ 1 - \f\ 2 - \Aa - \AI} 
where \ > \ 2 > A 3 > A 4 > are the eigenvalues of the 



matrix PabQAb witn QAB = ( a y ® °~y)P\ B {°~y ® ° 



'yh 



a Pauli matrix, and ® denotes the tensor product. The 
result of this calculation provides us with the expression 



C{Pab) 



2\S a \\Sf 



f aGb + Ga f b 



(23) 



This quantity is indeed restricted to the interval [0,1], 
as follows from the Cauchy-Schwartz inequality and the 
ine quality (^F ~G~ - ^F~G~f > 0, 2|5^||S S | < 
2y/ F A G A F B G B < F A G B + G A F B . The state described 
by Eq. (f2"2"| is trivial (i.e., not entangled or classical) only 
for zero overlap S A = and/or S B = 0. Physically, the 
vanishing of the overlap between the wave functions f(x) 
and g(x) in either of the two regions A and B implies, 
that these orbital states are perfectly distinguishable via 
a local measurement. In this situation the correspond- 
ing density matrix Pab ca n be written in a convex form 

PAB — J^iPiP^A ® P^S ' w ith probabilities pi > and 
y^ - Pi = 1, and thus is separable. On the other hand, for 
Sa b 7^ one cannot perfectly distinguish between the 
different orbital states via local measurements, result- 
ing in an interference between the different terms of the 
anti-symmetric wave function Eq. J21|) , a non-separable 
density matrix, and a finite concurrence. 

Next, we analyze the reduced density matrix in the 
outgoing leads, i.e., in the new basis {uu, ud, du, dd}. 
Exploiting the possibility of exchanging real space pro- 
jection and mixing, we can simply rotate the projected 
density matrix according to 



PAB 



= (17® U)p AB {U^ ®U r ). 



(24) 



The diagonal elements of the density matrix pab directly 
provide the detection probabilities P aia2 (AB), which 
then can be used in the calculation of the Bell correlation 
function Eab, Eq. (TJ|, 



Eai 



cos 2 {26) -sin 2 {26) C(pab) cos ipab, (25) 
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where the angle pab is given by the overlap integrals, 
<Pab = axg(S/i,Sg)] combining Eqs. ([25]) and (j2l)|) we 
immediately recover the original expression ([9]). Choos- 
ing four different intervals A, B, A', and B' (note that 
the selection of these intervals is non-trivial in the gen- 
eral situation discussed here, as p> ab now involves over- 
lap integrals), we can set up the Bell inequality §5§ 
and find the result expressed in terms of concurrences 
Ca'b — C(pa'b), 

| sin 2 (29) [Cab cos pab - Cab 1 cos p A B' (26) 
+Ca' b cos pa' b + Ca> b> cos p A > b'] + 2 cos 2 (29) | < 2. 

Choosing plane waves for f(x) and g(x), the concurrences 
take the value Cab = Cab 1 = Ca> b = Ca 1 B' — V and the 
Bell inequality reduces to the simpler form found earlier, 
see Eq. (ITU) ; for the general case, the degree of violation 
depends separately on the shapes f(x) and g(x) of the 
orbital wave functions in each region A, A' , B, and B' 
via the corresponding concurrences. 



D. Particles with spin 

So far, we have considered only spinless particles or, 
more exactly, two electrons in a spin-triplet state with 
the same spin polarization of the electrons, x+i ( a i > <?2) = 
^o-iT^T an d X-i( tJ i) (J 2) = $a 1 i5a 2 l- Since all spin- 
dependence of the Bell inequality is encoded in the over- 
lap Q of the spin wave- functions, see Eq. (0, one con- 
cludes that the above results are valid as well for the 
third maximally entangled triplet state, Xo( <7 i^ cr 2) = 
(<5o-iT<5o- 2 i + <5fTi i <5o- 2 i ) / \J~2 with Q = 1. On the other 
hand, the character of violation is modified for the spin- 
singlet state x sg (ci, c 2 ) = (8<rit&<ni ~ <W<*<73t)A/2 with 
Q = — 1. Choosing a set of optimal time intervals A, B, 
A' , and B' , the resulting Bell inequality takes the form 

£ = I V\/2sin 2 (2(9) - cos 2 (26»)| < 1. (27) 

The main difference to the previous result for spin-triplet 
states is that this inequality can be violated only for a 
sufficiently large visibility factor l/v2 < V < 1 and 
a beam splitter with a mixing angle 9 sufficiently close 
to optimal, 9 G [tt/4 — # c ,7r/4 + 9 C ], where the critical 
angle 9 C is given by sin 2 (2(9 c ) = 1/(1 + Vy/2), cf. Fig. M 
This result allows one to distinguish between triplet and 
singlet incoming states by measuring a Bell inequality 
involving only orbital degrees of freedom, see also Ref. @. 

Moreover, assuming that the incident electrons have 
opposite spin polarization, i.e., their spin state can be 
written as a superposition x in = ax® + /?X SS j the de- 
gree of violation of the orbital Bell inequality gives a 
lower bound on the value of the concurrence in the 
spin part of the wave function. Indeed, in this case 
Q(x m ) = \ a \ 2 — \P\ 2 and the maximal violation of the 
orbital Bell inequality for V — 1 and symmetric scat- 
tering is given by £ max = |Q|v2- At the same time, 




X 



FIG. 4: Spin-singlet source: The voltage pulse V(i) injects 
a singlet-pair of electrons into the lead s; the 7r/4 four- 
terminal splitter distributes the particles with equal proba- 
bilities among the two leads u' and a. The resonances in the 
quantum dots select the desired energies ei and £2; dots resid- 
ing in the Coulomb blockade regime inhibit the propagation 
of two electrons into the same lead, such that the two-particle 
incident state involves one particle in each of the leads u and 
d. 

the concurrence C(\ ln ) — |o; 2 — /3 2 1 > \Q\, with equal- 
ity established for real a and (3 (note that C(x in ) gives 
the degree of (useful) spin entanglement in the outgoing 
leads u and d). Hence, measuring the entanglement f ma x 
of the orbital part of the wave function (which leaves the 
spin component untouched), provides a (lower) estimate 
of the degree of spin entanglement of the incoming state. 
If the spin wave function of the incoming electrons fac- 
torizes, x m ( (J i> a i) — ^o-it^li the orbital Bell inequality 
never can be violated since in this situation the electrons 
are distinguishable and thus the detection of an electron 
with given spin in one of the outgoing leads always allows 
to determine its origin. 

III. SPIN-SINGLET /TRIPLET SOURCES 

Finally, we discuss the potential experimental real- 
ization of the proposed Bell test. The source of spin- 
entangled incoming particles can be realized with the 
help of a beam splitter, followed by leads with dots serv- 
ing as energy filters defined through resonance levels at 
energies e\ and £2, see Fig- [4] We assume that both dots 
reside in the strong Coulomb blockade regime; applying 
a single-electron voltage pulse to the source lead s, two 
electrons in a singlet state are detached from the Fermi 
se o 14 i 15 . There is only one scattering process, involving 
trajectories where the electrons tunnel through different 
quantum dots, for which the two electrons reach the sec- 
ond beam splitter. The incoming state in the leads u and 
d then is of a spin-singlet type with different energies E\ 
and £2 as defined through the dot resonances. All other 
scattering processes with only one or no electrons propa- 
gating towards the second beam splitter are irrelevant as 
they do not contribute to the correlation measurement. 

A spin entangled triplet state can be generated 
with the help of spin-polarized reservoirs with polar- 
izations | and I attached to the leads s and s, re- 
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spectively. Applying a single-electron voltage-pulse to 
each reservoir, two electrons with opposite spins are in- 
jected into the leads s and s, see Fig. 01 The state 

|*ss) = J dx!dx 2 f(x 1 )g(x 2 )%^(x 1 )i)l l (x 2 )\0) incident 
on the symmetric beam splitter emerges with a compo- 
nent 

|* a / a /) oc J dx ± dx 2 [g(x l )f(x 2 )tpl, l {x 1 )tpl, r [x 2 ) (28) 
+f(x 1 )g(x 2 )^(x 1 )$l, L (x 2 )]\0) 

describing electrons scattered into different leads u' and 
d'; it is this component which can propagate through the 
subsequent energy filter and contribute to the current 
correlators. The propagation of this component through 
the quantum dots results in an entangled spin-triplet 
state of the form given by Eq. (fTJ) with f(x) = exp(ifc]Xi) 
and g(x) — exp(ik 2 x). 

Above, we have considered an idealized situation where 
only two electrons are present in the system, while in a 
realistic situation one deals with electronic reservoirs at 
finite temperature. The associated equilibrium fluctua- 
tions then generate noise signals which are of the same or- 
der as the correlations associated with the injection of the 
two electrons. We note, however, that the correspond- 
ing equilibrium current correlators (I ai (x, ti)I a2 ( x i fe))eq 
assume significant values only for instantaneous or bal- 
listically retarded variables, i.e., at times t 2 = t\ and 
t 2 — ti = 2£/v F in the same leads and t 2 — t\ = 2£/v-p in 
opposite leads (here, I denotes the distance between the 
position of measurement and the reflecting dots). The 
Bell test involves correlations at time differences of the 
order of r = 2tt /5u> and a proper choice of the frequency 
mismatch Soj always allows one to render the contribu- 
tion from equilibrium fluctuations negligible. Another 
restriction on r is due to dephasing and electron-electron 
interactions; we then have to assume that the character- 
istic times associated with these processes are larger then 

T. 



IV. CONCLUSION 

We have discussed how to make use of quantum in- 
distinguishability as a resource to generate non-classical 
correlations: the indistinguishability of particles enforces 



proper symmetrization of their wave function and results 
in non-factorizable states. We have demonstrated how 
to generate such states with the help of quantum dots 
residing in the Coulomb blockade regime and have deter- 
mined their degree of entanglement as measured in a Bell 
inequality test based on auto- and cross-current correla- 
tors. In a real experiment, the latter are measured over 
a finite time or space domain. As a result, we obtain an 
interesting interplay between the measurement accuracy 
(time or space resolution) and the degree of non-locality 
as measured in the Bell inequality test: the more infor- 
mation is gained that locally distinguishes between the 
particles, the smaller is the degree of violation. Once 
the uncertainty principle allows for the identification of 
the particle, the Bell inequality cannot be violated any 
longer. This feature can be exploited in the design of 
experiments testing the above predictions: choosing a 
small energy difference Se — e 2 — £i allows for a slow 
measurement with a less stringent time resolution, while 
the violation of the Bell test remains observable. On the 
other hand, the energy difference 6e has to be chosen 
sufficiently large in order to avoid the influence of deco- 
herence or interactions. 

The above setup for spinless particles provides an alter- 
native for the observation of the two-particle interference 
as proposed by Samuelsson et aL— and recently observed 
by Neder et alZ; here, the role of the magnetic flux $ pen- 
etrating the Hanbury-Brown Twiss interferometer is re- 
placed by the time-delay of subsequent measurements in 
the correlator. Adding the spin degree of freedom, we are 
confronted with two distict situations: if the spin degree 
of freedom allows to distinguish between the particles 
(this is the case for the spin-state x{ cj It (j 2) = <5 CTl ^ <5 CT2 ) 
the Bell inequality is never violated. On the other hand, 
entangled spin states in the singlet or triplet sector (these 
are the states x sg an d Xo r ) can generate maximal viola- 
tion of the Bell inequality; finally, the superposition of 
these states reduces the spin-entanglement and the de- 
gree of violation in the orbital Bell inequality gives a 
lower bound on the spin-concurrence, with an ideal mea- 
surement providing the best bound. 

We thank Gordey Lesovik for discussions and acknowl- 
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dation, the RFBR grant No. 06-02-17086-a and the Pro- 
gramm "Quantum Macrophysics" of RAS. 



1 C. W. J. Beenakker, in Quantum Computers, Algorithms 
and Chaos, Proceedings of the Int. School of Physics "En- 
rico Fermi" Varenna 2005, Vol. 162 ( IPS Press, Amster- 
dam, 2006); |arXiv:cond-mat/0508488| 

2 N. M. Chtchelkatchev, G. Blatter, G. B. Lesovik, and Th. 
Martin, Phys. Rev. B 66, 161320(R) (2002). 

3 G. Burkard, and D. Loss, Phys. Rev. Lett. 91, 087903 
(2003). 

4 P. Samuelsson and M. Biittiker, Phys. Rev. B 73, 041305 



(2006). 

5 V. Giovannetti, D. Fmstaglia, F. Taddei, and R. Fazio, 
Phys. Rev. B 74, 115315 (2006); V. Giovannetti, D. 
Frustaglia, F. Taddei, and R. Fazio, Phys. Rev. B 75, 
241305(R) (2007). 

6 P. Samuelsson, E. V. Sukhorukov, and M. Buttiker, Phys. 
Rev. Lett. 92, 026805 (2004). 

7 I. Neder, M. Heiblum, D. Mahalu, and V. Umansky, Phys. 
Rev. Lett. 98, 036803 (2007) and I. Neder, N. Ofek, Y. 



8 



Chung, M. Heiblum, D. Mahalu, and V. Umansky, Nature 
448, 333 (2007). 

8 G. Burkard, D. Loss, and E. V. Sukhomkov, Phys. Rev. B 
61, R16303 (2000). 

9 A.V. Lebedev, G.B. Lesovik, and G. Blatter, Phys. Rev. 
B 71, 045306/1-9 (2005); note that with only two incident 
electrons there is no restriction on the accumulation time, 
see Sec. Ill B. 

J.F. Clauser and M.A. Home, Phys. Rev. D 10, 526 (1974). 

1 A. Aspect, P. Grangier, and G. Roger, Phys. Rev. Lett. 



49, 91 (1982). 

G. B. Lesovik, JETP Lett. 49, 592 (1989); M. Biittiker, 

Phys. Rev. Lett. 65 2901 (1990); Ya. M. Blanter and M. 

Biittiker, Phys. Rep. 336, 1 (2000). 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

L. S. Levitov, H. W. Lee, and G. B. Lesovik, J. Math. 

Phys. 37, 4845 (1996). 

J. Keeling, I. Klich, and L. S. Levitov, Phys. Rev. Lett. 
97, 116403 (2006). 



